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Summary 

The Langley Research Center has designed a 
swept, supercritical airfoil incorporating laminar- flow 
control (LFC) for testing at transonic speeds. Ana- 
lytical expressions have been developed and an eval- 
uation has been made of the experimental section 
drag, composed of suction drag and wake drag, by 
using theoretical design information and experimen- 
tal data. 

The analysis shows that, although the sweep- 
induced boundary-layer cross-flow influence on the 
wake drag is too large to be ignored and there is 
not a practical method for evaluating these cross- 
flow effects on the experimental wake data, the con- 
ventional unswept two-dimensional wake drag com- 
putation used in the reduction of the experimental 
data is at worst 10 percent too high. 

Introduction 

The total drag of laminar-flow-control airfoils 
may be determined by combining the measured suc- 
tion power (suction drag), required to maintain a 
laminar boundary layer, with the measured stag- 
nation pressure losses in the wake just downstream 
of the trailing edge (wake drag). The combination 
of the two drag forces (wake drag + suction drag) 
thus represents both the resistance of the air through 
which the airfoil moves and the penalty incurred in 
creating a suction flow through the airfoil surface. 
This total drag must be less than the wake drag of a 
comparable nonsuction airfoil in order to justify the 
use of boundary-layer suction. 

The suction drag is derived from considerations 
relating to the power required to Create suction 
through the airfoil surface. It is not an actual phys- 
ical drag acting to oppose the motion of the airfoil 
through the air, but a drag computed from suction 
power requirements. Not considered in the analysis 
presented in this paper are the penalties due to the 
weight and maintenance of the suction apparatus and 
the inefficiency of the suction compressor. 

The axisymmetric suction nozzle used in these 
tests was developed by Dr. Werner Pfenninger dur- 
ing early laminar-flow-control research. The profile 
shape of the nozzle and the calibration curve for this 
nozzle, given in the appendix, were communicated in- 
formally by Dr. Pfenninger during the design phase 
of these investigations. The original derivation of the 
equation for suction drag (eq. (13)) is also due to 
Dr. Pfenninger. 

The wake drag may be determined either by the 
integration of momentum deficits across the wake as 
computed from measured wake-rake pressures in the 


conventional manner of references 1 and 2, for two- 
dimensional airfoils, or by more sophisticated meth- 
ods that account for the airfoil sweep by analysis 
of the airfoil theoretical boundary layer. When the 
trailing edge is swept, the conventional wake-rake 
drag computation can be in error because of either 
inviscid cross-flow (defined as the deviation of the 
boundary-layer edge flow from the streamwise direc- 
tion) effects or viscous cross-flow (deviation of flow 
in the boundary layer from the edge direction) effects 
(refs. 3 to 6). Therefore, an error analysis of the wake 
drag is required for swept airfoils, since the mass-flow 
rate on which the momentum deficit is based can be 
affected by boundary-layer cross flow and resulting 
pressure gradients. 

The object of this report is the development and 
application of analytical methods for the evaluation 
of the coefficients of suction and wake drag for the 
slotted and porous laminar-flow-control (LFC) air- 
foils that have been tested in the Langley 8-Foot 
Transonic Pressure Tunnel (8-ft TPT) at transonic 
speeds over a range of chord Reynolds numbers. An 
overview of these experiments is presented in refer- 
ence 7, and basic experimental results for the slotted 
configuration are presented in reference 8. 

Even though the conventional wake-rake drag 
computation method does not take into account 
cross-flow effects, it was the method used to com- 
pute wake drag from experimental measurements in 
the LFC experiment of reference 8. Drag computa- 
tions based on the theoretical boundary layer on a 
swept airfoil are evaluated and presented as an indi- 
cation of the magnitude of the cross-flow effect. 

Symbols 

b airfoil span, measured perpendicular 

to tunnel centerline 

bjy airfoil span, measured along airfoil 

sweep line 

C p local pressure coefficient, 

(p ~~ Poo)/ Qoo 

Cq nondimensional local suction 

coefficient, -(pw) W8 /(pU ) oo 

c airfoil chord, measured parallel to 

tunnel centerline 

c d section drag coefficient 

c d s suction drag coefficient 

c d,w wake drag coefficient 

Q, total t° ta l drag coefficient, c ds + c d yy 


C\ 

section lift coefficient 

airfoil chord, measured normal to 
airfoil leading edge 

w 

velocity component perpendic- 
ular to free-stream direction (in 

CjV 


Z-direction) 

Cp 

specific heat at constant pressure 
specific heat at constant volume 

w 

a N 

velocity component in z-direction 
suction-nozzle flow coefficient 

Cy 

D 

drag 

7 

ratio of specific heats, Cp/c v (1.4 for 
air) 

D s 

suction drag 

A 

incremental amount 

d 

diameter 

6* 

boundary- layer displacement 

F 

area of duct in plane perpendicular 


thickness 


to suction nozzle axis 

e 

boundary- layer momentum 

f 

cross-sectional area of nozzle throat, 


thickness 


nd 2 N /A 

A 

sweep angle 

H 

boundary-layer shape factor, 8* /0 

P 

density 

L s 

suction power, ^jr 

See figures 2 and 4 for geometric orientation of 

M 

Mach number 

following coordinates: 

m 

mass, slugs 

X 

coordinate in flow direction parallel 
to tunnel centerline 

rh 

mass flow, slugs /sec 

X 

coordinate perpendicular to leading 

P 

pressure, psf 


edge 

Q 

dynamic pressure, psf 

Y 

coordinate orthogonal to X in 

R 

gas constant, c p — c v 


airfoil reference plane, positive up 

Rc 

free-stream Reynolds number based 

y 

coordinate parallel to leading edge, 


on chord 


positive up 

Rd 

Reynolds number based on suction- 

z 

coordinate orthogonal to X and Y 


nozzle throat diameter 

z 

coordinate orthogonal to x and y 

s 

airfoil surface area, be 


(also used for local surface normal, 
positive out from surface) 

s 

distance along surface 

Subscripts: 


T 

temperature, °R 

e 

potential flow at boundary- layer 

^Tt,sc 

temperature rise of suction process 
(eq. (5)) 

i 

edge 
ith term 

U 

velocity component in free-stream 
direction (in X-direction) 

N 

nozzle or normal 

u 

velocity component normal to 
leading edge (in ^-direction), 
U cos A 

P 

boundary-layer velocities parallel to 
free-stream direction at boundary- 
layer edge 

V 

velocity component parallel to 

s 

suction 


leading edge (in y-direction), 
U sin A 

sc 

suction chamber 
stagnation condition 


t 


total combined suction and wake effect 


2 



tr boundary- layer transition 

u , v boundary-layer thickness definition 

velocities (eqs. (19b) and (32)) 

W wake 

W.R. evaluation at the wake-rake 

position, X/c = 1.094 

ws airfoil surface 

oo free stream 

Abbreviations: 

L.E. leading edge 

T.E. trailing edge 

LFC laminar-flow control 

TPT Transonic Pressure Tunnel 

2-D two-dimensional 

Experimental Apparatus and Model 

A sketch of the LFC airfoil configuration, which 
incorporates laminar-flow control into an advanced 
supercritical airfoil, and the shock-free design pres- 
sure distribution for this profile are shown in fig- 
ure 1. The three types of boundary-layer instabil- 
ities that were taken into account during the design 
process are noted on the figure. This airfoil has 
a nearly full-chord discrete suction surface with an 
internal suction-flow ducting system. The continu- 
ous surface suction as modeled in the design theory 
can be approached only asymptotically by the ac- 
tual suction surface of the wind-tunnel model. For 
the LFC experiment in the 8- ft TPT (ref. 7), two suc- 
tion surfaces were evaluated: one with narrow (about 
0.003 in.), closely spaced spanwise slots and one with 
spanwise strips of perforated titanium about 0.5 in. 
wide and 0.4 in. apart (see fig. 2(c)). The slotted 
surface suction concepts are discussed in more detail 
in references 7, 8, and 9, and the porous surface suc- 
tion concept is discussed in reference 10. Sketches of 
the design laminar regions for these LFC airfoils are 
shown in figure 2. The laminar regions are bounded 
by turbulent wedges that develop on the airfoil; these 
wedges are generated by the junctures of the tunnel 
liner wall and the airfoil. Suction extends rearward 
to X/c = 0.96 on the slotted- upper-surface airfoil 
and to X/c = 0.89 on the porous-upper-surface air- 
foil. On the lower surface, only slotted suction was 
used; the suction extends from near the leading edge 
to X/c = 0.842. The model was mounted vertically 
in the wind tunnel with 23° of sweep and spanned 
the test section from ceiling to floor. 


Figure 3 shows the chordwise suction distribu- 
tions theoretically determined to maintain full-chord 
laminar flow over a range of Reynolds numbers for 
the upper and lower airfoil surfaces. The theoreti- 
cal suction distribution for the design case (M 0 0 = 
0.82, R c = 20 x 10 6 ) — derived from a linear, paral- 
lel, incompressible boundary-layer stability analysis 
(refs. 11, 12, and 13) — is also presented in table I 
along with the corresponding suction drag. The the- 
oretical chordwise pressure distribution (fig. 1) was 
used as input to boundary-layer stability codes to de- 
termine these suction requirements. The expressions 
used to compute the suction drag attributable to ei- 
ther a theoretical or an experimental suction distri- 
bution are derived in a following section. 

Discussion of Analytical Methods 

In order to analyze suction and wake drag, the 
swept LFC airfoil is considered to have constant cross 
section and infinite span in steady transonic flow. 
The local suction velocity through the surface is as- 
sumed to be constant along sweep lines parallel to the 
leading edge. The flow field over the airfoil is sub- 
divided into a “suction flow” containing stream- 
lines entering the airfoil suction surface, a “wake 
flow” containing streamlines entering the trailing- 
edge wake, and a “potential flow” containing all other 
streamlines. The suction system (fig. 2) consists of 
surface slots or perforations, internal airflow meter- 
ing, ducting, and the compressor. The suction flow 
passes through the suction system and is returned 
to the wind-tunnel circuit far enough downstream so 
as not to affect the airfoil wake. An “energy de- 
fect” or decrease in total pressure occurs as a result 
of viscous dissipation in the sucked portion of the 
boundary layer and because of flow through the air- 
foil surface into the airfoil ducts. The portion of this 
energy defect that occurs upstream of the suction 
nozzle throat is represented as the suction drag, c<i s . 
Total drag is then defined as the sum of the drag as- 
sociated with the energy defect in the suction system 
and the wake drag associated with momentum defect 
in the wake: 

c d, total = c d,s + c d,W 

Suction Drag Coefficient 

The expression for the suction drag coefficient de- 
rived below is that used throughout the LFC tests 
in the 8-ft TPT. It is important to emphasize that 
this suction drag represents the energy deficit due 
to moving air through the airfoil surface and inter- 
nal airfoil ducting as far as the throat of the suc- 
tion nozzle (fig. 2(c)) and then returning this air to 
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free-stream temperature and pressure. No attempt 
is made to evaluate energy losses in the connector 
hoses, suction controls, or compressor. 

Derivation . The net power delivered to the 
suction flow, L 5 , is defined as the product of the free- 
stream velocity and an equivalent suction drag: 

L$ — DgUoo 

The suction drag coefficient is then 

D $ D s 2 D s 

Qoc S 0. 5 poo S poo S 
__ ‘ZDsUqq _ 2 L s 

Note that this equation applies to one surface 
(area = S) of the airfoil. The suction power L s 
equals the energy defect per unit time resulting from 
the flow through the internal airfoil airflow ducting 
system. The suction power L s = c p mATt^c (where 
c p is the specific heat at constant pressure) has pre- 
viously been defined in reference 14 (ch. 5, eq. (5.7)). 
Therefore, since fn = -^ = [— pwS\ ws , where the 
negative sign arises from the definition of z , the local 
surface normal coordinate, as positive out from the 
surface 


pressure is given (in the isentropic perfect gas 
approximation) by 


A T< 


compression 


Tt.s 



( 6 ) 


In addition to being compressed to free-stream pres- 
sure poo, the air (essentially at rest in the duct) must 
be accelerated to free-stream velocity, [Too* Using the 
adiabatic energy equation, this further temperature 
increase is 


^^acceleration — 


7~1 

2 


M^Too 


( 7 ) 


Upon substitution of equations (6) and (7) into equa- 
tion (5) 




2 1 m 2 t 



( 8 ) 


Substitution of equation (8) into equation (4b) for 
the suction drag gives 


L s — f A Tt, sc CpS[ (pw)u;s] d ^ ^ ^ — DsUqq (2) 

J L.E. 

or 

/■T- E • 

L s = I ATt t)S c c p^Poo u oo[~ {p w )ws/ Poo u oo] ( “ ) (3) 

•/L.E. 

Now, by substitution of equation (3) into equa- 
tion (1) the suction drag coefficient may be expressed 
as 

2 Jle ATf jSC CpS'poo w oo[ — {p w )ws/ Poo u oo] d (c ) f A ^ 
Cd ’ s = ^ ( } 



It can be shown that since, in the standard perfect 
gas approximation, the free-stream speed of sound 
Uoc = y/^RToc- and 7 = c p /c v , R = c p - c v , and 
Moo — Uoo! a oo 


or 


T E 

Cd,s = 2 < ^3°° J AT t<ac [-{pw ) ws / P00U00} d(^) ( 4 b) 

The term A T^ sc is the temperature increase involved 
in returning air in the suction duct to free-stream 
pressure and velocity. 

m ,sc — A ^compression + AT acce i era tj on (5) 

The temperature increase involved in compressing 
the air from suction duct pressure p sc to free-stream 


U*o = 2 ( 10 ) 

Thus with Uqo = Uoo cos A we can substitute equa- 
tion (10) into equation (9) and write the term outside 
the integral in equation (9) as 


‘ICpUoo 

C/3 “ 

u oc 


2c p Uoc cos A 

2 

2 cos A 

M2,( 7 -l)T 00 


( 11 ) 
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Equation (9) becomes, therefore, 



or, finally, the equation used to evaluate the suction 
drag coefficient from experimental data 


_ cos A /" TE - ^ 

° d ' s ~ OlMg Jl.E. Q ' SC 


X 



where Cq, sc = 




The evaluation of the term 


Cq,sc “ 


' ~{pw)ws 
. Poo u oo . 


from experimentally measured quantities is presented 
in the appendix. Note that, for values of sjc where 
there is no suction, Cq sc is zero and so the integrand 
in equation (13) is also zero. 

Equation (13) allows the evaluation of the equiv- 
alent suction drag coefficient for one surface of the 
LFC airfoil due to the energy deficit associated with 
returning suction duct air to free-stream conditions. 
Although the energy lost in forcing air through the 
skin slots or perforations is accounted for, the energy 
dissipated in the wall boundary layer of the suction 
hoses downstream of the nozzle and the efficiency of 
the suction compressor are not represented in any 
way in this equation. 

Note that what is computed in equation (13) is 
a section drag coefficient, since the integration is 
carried out only in the stream direction, and the 
flow through the surface takes place over an area 
(chord x unit span) of model surface. As the wake 
drag coefficient computed from a wake rake is also 
a section drag coefficient, the drag represented by 
equation (13) is the appropriate suction drag to 
add to the wake drag to obtain the total drag of 
the suction airfoil for comparison with the drag of 
nonsuction airfoils. 


Sample computation . As an example of the use 
of equation (13) in the data reduction process for the 
slotted airfoil used in the LFC experiment, the terms 
of the equation are computed for duct 9 (a typical 


duct in the center of the upper surface) for a typical 
case at design Reynolds number (R c = 20 x 10 6 ). 
Taking the terms in order from left to right 

A = 23° 

Moo = 0.8188 

T*,oo - 547.54°R 

Poo — 937.3 psf 

and thus (from the standard adiabatic perfect gas 
equation for M 0 0 and Tt )00 ), T 0 0 = 482.80°R. 

The limits of integration are, for this example, the 
upstream and downstream boundaries of the duct: 
A (s/c) = (30.969 in. - 25.571 in.)/78.096 in. = 
0.06912 and the quantities in the terms under the 
integral are assumed constant over the chordwise 
extent of the duct. 

For measured values of nozzle throat drop A p^ — 
—29.28 psf, duct temperature Tt )SC = 532.17°R, and 
duct static pressure p sc = 433.98 psf, Cq sc for 
duct 9 may be calculated from equation (A3) of the 
appendix to be 


Cq^sc “ 


' -(pw)ws ' 
. PoqUqc . 


0.000122285 


With 7 = 1.4, we may now compute (from 
eq. (13)) that the contribution of duct 9 to the suc- 
tion drag coefficient is 0.00002555. Although equa- 
tion (13) is derived as an integral, in actual use it is a 
discrete summation, and this suction drag coefficient 
computed for duct 9 is the 9th term of the summa- 
tion. There would be 24 parallel terms for the slotted 
airfoil upper surface (21 terms for the perforated) and 
another 21 terms for the lower surface. 

It should be noted that the numbers in this ex- 
ample are given to greater precision than is actually 
available from the instrumentation. In general, pres- 
sure errors are at least on the order of 1 psf (except 
for the instruments used to measure tunnel total and 
static pressures, which are accurate to 0.2 psf, giving 
a Mach number precision of 0.0007) and temperature 
errors at least 1°R. Based on these estimates, the 
number computed as Cq, sc of duct 9 (0.000122285) 
has a possible error of ±0.000002 {Cq, sc = 0.000120 
to 0.000124), and the computed suction drag co- 
efficient has a possible error of ±0.000001 (duct 9 
A c d , s = 0.000025 to 0.000026). 


Wake Drag Coefficient 

Five expressions for the wake drag of a swept 
2-D airfoil are derived, three based on theoretical 
boundary-layer parameters that would be difficult to 
measure experimentally, and two based on wake-rake 
data that do not fully account for the cross-flow ef- 
fect induced by the airfoil sweep. Note that the three 
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methods (the “leading-edge-perpendicular” method, 
the “streamwise” method, and the “semiempirical” 
method) based on theoretical boundary-layer param- 
eters assume that these boundary-layer parameters 
can be carried across the trailing edge singularity into 
the wake. 

Baals and Mourhess method . The wake drag 
of the airfoil in the LFC wind-tunnel experiment 
(ref. 7) is evaluated by the conventional momentum 
deficit method of reference 2. This is the method 
that is commonly used for the computation of the 
wake drag of unswept airfoils. The design calcu- 
lations for the slotted LFC airfoil in the 8-ft TPT 
(ref. 15) indicated that the wake outer edge cross- 
flow angles at the wake rake are within the angle 
range for total pressure measurement error less than 
1 percent, that is, ±20° with respect to the wake-rake 
total tube plane. The wake rake is located 8 in. or 
9.4 percent chord downstream of the airfoil trailing 
edge and aligned, in the airfoil-chord plane, parallel 
to the tunnel reference centerline. The rake face is 
aligned (in the downflow plane) perpendicular to the 
camberline tangent plane through the trailing edge of 
the airfoil with 0° flap deflection. (See figs. 4 and 5.) 

Flow visualization data, using tufts on the rake, 
indicate that, while the extreme end of the rake above 
the airfoil upper surface is aligned with the flow and 
the extreme end of the rake below the airfoil lower 
surface is at about -10° yaw (that is, in the airfoil 
chord plane) to the flow, the center of the rake may be 
inclined as much as 4-30° yaw to the flow. This will 
result in a computed wake drag that is 5 to 10 percent 
high (1/2 to 1 count of drag at design Mach number 
and chord Reynolds number of 10 x 10 6 ) based on the 
total-tube angle error and wake drag calculations of 
references 16 and 17. These results indicate that it 
would have been extremely difficult to devise a better 
wake rake, since the local spanwise flow angle seems 
to vary by 40° over the thickness of the wake at this 
position downstream of the swept airfoil. 


Leading- edge-perpendicular method. Realiz- 
ing the limitations involved in application of the stan- 
dard 2-D wake drag computation to a swept airfoil, 
we will present here an analysis that will give some 
indication of the relative magnitude of the cross-flow 
effect on the wake drag. 

Applying the approach used by Raetz (appendix I 
of ref. 5), we may begin with the following expression 
for the wake drag of an untapered swept airfoil. 
The wake drag is defined as the momentum loss of 
the external airflow far downstream of the trailing 
edge where the static pressure equals the free-stream 


value: 

2 D\y 

Cd ' W = <14o) 

or, in terms of the wake momentum thickness (see 
ref. 1), 

26 , v 

c d,W - — (14b) 

Now consider a spanwise segment of the airfoil, 
dy, with chord c/y. The velocity deficit for this seg- 
ment is determined far downstream of the trailing 
edge so that the static-pressure field everywhere in 
the wake is essentially equal to free-stream static 
pressure. Considering the sketch in figure 4, an ex- 
pression for the compressible flow per unit time of the 
momentum component in the streamwise direction, 
pU = p(u cos A + usinA), (through a unit area face 
of control volume just downstream of and parallel to 
the trailing edge) is 

^downstream = p{u COS A + v sin A )u dy dz (15) 


and the momentum flow 

PooUoo = pooi^oo cos A + Vqq sin A) 

in through the upstream face of the same control 
volume is 

^upstream = Pooi^oo cos A 4- t>oo sin A)Uqo dy dz (16) 

Thus, the integral of the difference between expres- 
sions (16) and (15) over the upstream and down- 
stream faces of the control volume equals the net flow 
per unit time into the control volume of the momen- 
tum component in the streamwise direction, which is 
the drag: 


/ oc r oo 

/ (^upstream ~ ^downstream) 

-oo J — 00 


With conservation of mass in the control volume, 


171 control volume — 



(pu - PooUoo) dy dz = 0 


Thus m contro i volume u oo COS A = 0 and ^control volume 
t^oo sin A = 0. Taking the sum of these zero terms 
with the difference of expressions (16) and (15) 


D w 


/ OO roc 

/ (^upstream ~ ^downstream) 

-oo J — OO 


4" ^control volume ^ oo COS A 
4“ ^control volume ^oo sin A 
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and rearranging the terms, the drag can be expressed are 


D w 



(u 0 o — u) cos A 


{@uu) c 


-f. 


pu 


oo Poo u 00 



dz 


+ (Voo - v) sin A 

dy dz 

(g ) - f 00 p u ( 

!_jA 



U °° J-O 0 Poouoo V 

Voo) 


> 



(19a) 


Since all variables are assumed constant along sweep 

lines, f dy = btf and, because we want to compute a far downstream of the trailing edge. Equation (18b) 
section drag, bj^ = 1 and so may be expressed as 


/oo 

pu[(uoo — u) cos A 

-oo 

+ (^oo — v) sin A] dz (17) 


or 


D w 


O.bpooUlc cos A 

2 r°° 


J 

J - 0 


PooUlc J_ 0C 

2>Poc u oo 
PoqUoqC 


pu (u 0 o — u) + (^oo — v ) tan A 


dz 


Voq tan A 
Uoo 


u oo f 

Uoo J. 

f. 


pu 

Poo u oc 


(1 — u/uoo) dz 


pu 

oo Poo u oo 


(1 - v/voo )dz 


c d,W _ 2uo 
cos A Ur 


00 


Up o( g uu) c 
UooC 


. V 00 \Vvu ) oo , A 

+ — tan A 


Uc 


00 


(19b) 


Letting Uoo = Uoo cos A and Voo = Uoo sin A, 
substitution of equations (19a) into (19b) gives 


c d,W 
cos A 


= 2 


+ 


{@uu) 


C 

(0vu) 


— cos 2 A 


— sin A cos A tan A 


or 


For S = 6c, the wake drag coefficient (eq. (14a)) may 
be expressed as 


Cd,w = 2 cos A 


\yuu)oc 2 a 
— cos A 


2 D w 

W 


(18a) 


+ 


^ksin 2 A 

c 


( 20 ) 


and so 


c d,W _ 2/300^00 
cos A PooUqqC 


+ 


Uoo [* 

UooJ-o 

j: 


pu 


oo J— oo Poouoo 
Voo tan A f°° pu 


(1 — u/uoo ) dz 


Uoo J— oo Poouoo 
X (1 — v/voo) dz 

The components of the momentum thickness 

e m =r — ( 1 --) dz 

J- oo Pe u e V u e ) 


(18b) 


and 


= f°° _py_ ( 
J — oo Pe^e V 


1 ] dz 

v e > 


which is the wake drag coefficient for a constant- 
chord swept airfoil according to Raetz (appendix I 
of ref. 5). In essence, A represents the sweep angle 
of the isobars and is used to determine the direction 
of the component of velocity that is modified by the 
pressure gradients downstream of the airfoil trailing 
edge. 

The momentum thickness terms in equation (20), 
(0 UU ) oo an d (0vu) o o’ are momentum deficits far down- 
stream of the airfoil in the direction perpendicular to 
the leading edge. To evaluate (O uu ) 0 0 and (6 vu )o o in 
terms of quantities at the airfoil trailing edge, one 
must consider the momentum thickness of the wake 
at the trailing edge with cross-flow effects due to 
sweep. Again following the approach used by Raetz 
(appendix I of ref. 5), and by Squire and Young 
(ref. 1), for the momentum integral equations down- 
stream of the trailing edge (x- and z-directions with 
v constant) and assuming that the shearing stress is 
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zero, 


+ ( tf + 2 )^i^ = 0 (21a) 

Wg dx 

where H = 6*/6. Because the net transport of 
spanwise momentum into a control volume aligned 
perpendicular to the trailing edge must be zero from 
the assumption that all parameters are invariant 
along sweep lines, 

— (u e 0t, u ) = 0 

and, therefore, 

u e 0 vu = Constant (21b) 


H). Both In {uoo/ue) and H simultaneously decrease 
monotonically along the wake from the trailing edge 
to a large distance downstream, i.e., 

ln(uoc/u e ) = 0 (far downstream, where H — 1) 

If we assume, as in reference 5, that 

_ Constant 

H — 1 


then the integral in equation (25) can be written 



where x is the perpendicular distance from the trail- 
ing edge and u e is the velocity (perpendicular to the 
trailing edge) at the edge of the wake. Equation (21a) 
may then be written (dividing by 

@uu) 


@uu dx dx \ Ut(X) J 


(22) 


Integrating equation (22) over the distance x from 
the trailing edge to a point sufficiently downstream 
such that u = u 0 0 and p = Poo 


In 0-nu 

T F 

' = - {H + 2) In — 


00 Woo 


T.E. 


+ / TE ' ln (23) 

J oo u oo dx 


Substitution of equation (26) into equation (25) gives 



Furthermore, from equation (21b), j^{u e 9 vu ) = 0, 
and so we can integrate exactly to get 


(6U)qq _ »e,T.E. 
(^uu)x.E. u oo 


(28) 


For large distances downstream of the trailing edge, 
H = 1 at x = oo, so 


In 


(^uu)t.E. 

(@uu) oo 


+ (#T.E. + 2) In 

u oc 

t T.E. 

= / In — — dH 

JH= 1 Wqo 


(24) 


i?X.E. IS the value of the shape factor at the trail- 
ing edge and may be obtained from boundary-layer 
calculations. Then in terms of the integral in equa- 
tion (24) and u e ^. e. 


Then, by substituting equations (27) and (28) into 
equation (20), we get the expression for the coeffi- 
cient of wake drag 

p +5 

2cosA[ /«e,T.E.\ 2 ' 2 

c d,W ~ I (^uu)t.E. ^ — J COS A 

+ (^vu)t.e. (““ ) sin 2 a| (29) 

Equation (29) allows the computation of a wake 
drag that can be added to the suction drag of equa- 
tion (13) to yield the total drag of the LFC airfoil: 


(Moo _ ^ e ,T : E_ 1 \ J?TE+2 e / 1 TE ln^dH (25) 

(duu)x.E. V u °° ' 

Evaluation of the integral in equation (25) re- 
quires a relationship between the static pressure in 
the wake (determines the value of u e ) and the veloc- 
ity distribution in the wake (determines the value of 


c d, total = c d,s + c d,W ( 3 °) 

To illustrate the use of equation (29), theoreti- 
cal values of the boundary-layer momentum thick- 
nesses and edge velocities computed for the LFC air- 
foil at design conditions (taken from unpublished re- 
sults based on the methods of ref. 11) presented in 
table 11(a) are used to calculate and the results 
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are given in table 11(b). Note that the boundary- 
layer parameters 0 and 8 * based on flow perpendicu- 
lar (uu) and parallel (vu) to the leading edge of the 
swept airfoil are of roughly the same order of mag- 
nitude and that the second term in equation (29) is 
as large as 22 percent of the first term. In the origi- 
nal theoretical calculations for the design of the LFC 
airfoil, the onset flow was broken down into orthog- 
onal components perpendicular and parallel to the 
airfoil, so that equation (29) is the appropriate com- 
putation of the drag from the theoretical boundary 
layer. This equation has not been used for the ex- 
perimental drag coefficient, however, because of the 
difficulty in measuring velocity components u and v. 


so that 



Or, dividing equations (33) and (34) by 6 


Streamwise method . The computation rep- 
resented by equation (29) can be carried out with 
theoretical values of 0 UU and 6 VU obtained from a 
boundary-layer calculation in a plane perpendicular 
to the airfoil leading edge. However, such boundary- 
layer values cannot be measured experimentally. 

In an effort to evaluate the cross-flow effect on 
wake drag from quantities in the streamwise direc- 
tion, which could be measured experimentally, the 
following approach due to Groth (ref. 3) was used. 
New coordinates in the directions parallel and per- 
pendicular to the free-stream flow are defined and 
the local boundary-layer velocity coordinates u and 
v transformed into u p and v p as follows. (See fig. 4.) 



Thus, after substitution of equations (35a) and (35b), 
the bracketed terms in equation (29) may be ex- 
pressed as 


u e 


v e 


Uv = — U + jf'V 
' e 


P U e " Ue 


(31) 


V e 


He 


Vp = u+ —V 


U< 


Ue 


Note that this analysis applies to only one side of the 
airfoil at a time. 

The boundary-layer displacement and momen- 
tum thicknesses are defined in the new coordinate 
system as follows: 



OiipUp — 

P VpUp 

&VpV p = 


r (l _ PU P 
Jo \ PooU e 

Jo PooUe 

r-eat-fi- 

J 0 PooU e \ 

f°° P V P U P dz 

Jo Poouy z 

L 


00 PVp 
PooUi 


dz 




dz k 


(32) 



where all the terms on the right-hand side of the 
equation are evaluated at the trailing edge (T.E.). 

If we now assume (see fig. 4) that the velocity 
along the sweep line at the edge of the boundary 
layer at the trailing edge, (v e ) T E , can be taken as 
equal to the component of free-stream velocity in this 
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direction, Vqo, then 


-'OC t^oo 


u 0 o tan A 


\ u e / t.E. u e,T.E. ^e,T.E. 
and so the right side of equation (36) can be written 


^ UpUp "h 


(@v p Up ~ Kp') t an A 


x (jfL > \ a?5 C0S 2 A+ f M sin 2 A 


+ ^ ®vpv p tan A 

r / x ^+3 

f U e \ Z 

x sin A cos A 

\UooJ 


sin A cos A 


„ (UeoXV 2 

V, J “» A 


[(^^pup — fi Vp ) tan a| 

/ U» \ ^ o 


u e\~ 2a 

— cos A 


®V p Vp 


sin 2 A 


sin A cos A 


+ (OupUp + 0 V p V ^j (^) Sln2 A 

+ (dvpUp - Kj) tan A sin 2 A 

We can now substitute these expressions into equa- 
tion (29) to obtain (arranged as terms in 0 UpUp , 9 VpUp , 

^t? n 5 an d 


2 cos A f / u e 2 » 

W = T~ (— ) cos A 


sin 2 A 



u e 

, u oo 

+ ®v p u p 

-( 

u e 

V 

. U QO 


— — ) sin A cos A + tan A sin 2 A 


6* ( — I cos 2 A + sin 2 A tan A 


sin * A > (37) 

U OG 


The boundary-layer parameter components (terms in 
9 and 5*) of the wake drag in this form may now be 
evaluated with theoretical values from unpublished 
results using the methods of reference 12 as presented 
in table 111(a). These are the same basic theoretical 
values as were used in the leading-edge-perpendicular 
method, transformed into the streamwise coordinate 
system. It can now be seen in table 111(b) that in 
the streamwise analysis, the contribution of spanwise 
terms at the airfoil trailing edge may be considered 
small, as the largest cross-flow term (the third term 
in eq. (37), a multiple of 8* p ) is 0.2 percent of the 
first term. 


Semiempirical approximation . In the semi- 
empirical approach for swept wings of reference 1, 
the wake drag coefficient is expressed in terms of 
variables at any point x downstream of the trailing 
edge as 

H x + 5 

20 x ( U x \ 

^ - ~ k) 

Taking this point x as the position of a wake rake, 
we may write 


Hw+ 5 

2 0\y ( U e ,W \ 2 


c d,W = 


The momentum and displacement thicknesses as a 
function of the ratio of the measured resultant ve- 
locity in the wake to that in the free stream may be 
determined from the expressions 


e w = J_ 

S ” = L 


e _pU_ 

-e PeV ti 


where U is the velocity determined at a wake rake 
oriented in the streamwise direction. The rake size 
and the downstream distance from the trailing edge 
are designed so that there is a point on each end of the 
rake such that U = £/ e , beyond which the integrand 
is zero. Distributions of both static pressure p and 
total pressure pt are measured on the rake and, using 
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standard compressible flow relationships, 


U 2 


27 P 
7-1 P 



(40) 


The velocity ratio may be expressed as 


U_ 

U e 


7-1 

P[Pe (Pt/P) 7 ~ 1 

P,Pe ( Pt,e/Pe ~ 1 


(41) 


where the density ratio in terms of the velocity ratio 
at the rake is 


p_ = 1 + ^M 2 
Pe ~ 1 + 2^1 Af| 

In terms of Mach number, 


This equation provides the velocity ratio that could 
actually be measured by a suitable rake or survey 
probe at the trailing edge in terms of the perpen- 
dicular edge velocity ratio. For the purpose of com- 
paring equations (29) and (38), equation (43) has 
been used to compute a streamwise velocity ratio 
from the theoretical perpendicular velocity ra- 
tio U ^ E ' and the result is presented in table II. 

The computation of wake drag from equation (38) 
using the theoretical boundary- layer values of ta- 
ble II for comparison with the value obtained from 
equation (29) assumes that the momentum thickness 
and shape factor (which are boundary-layer values) 
of equation (29) are adequate to this order of ap- 
proximation for use in equation (38) (a wake equa- 
tion). Note also that, for the purposes of this com- 
parison, the wake drag coefficient computation from 
equation (38) must be carried out separately for the 
upper and lower surfaces. 


U _ M 1 + 

U e ~ M e Aj 1 + V M2 


(42) 


The coefficient of equation (38) can now be 
computed from pressure data measured by a rake 
in the airfoil wake. This has been done for the 
slotted airfoil data at chord Reynolds numbers of 
10-15 million. A sample pressure distribution and a 
computed velocity distribution on the wake rake are 
presented in figures 6 and 7. The drag coefficients 
are presented in figure 8. For comparison, we also 
present in this figure the corresponding 2-D drag 
coefficient computed from the same wake-rake data 
by the conventional Baals and Mourhess method of 
reference 2. 

To compare the wake drag coefficient obtained 
with the semiempirical approximation of equa- 
tion (38) with that obtained by the leading-edge- 
perpendicular method of equation (29) or the stream- 
wise method of equation (37), using the streamwise 
theoretical boundary- layer values of table III, assume 
the inviscid relationship ^ e ,T.E. — Uoc sin A. Then 


w e,T 




E yju j>.E. (Uoo sin A) 2 
Uoo cos A 

(lf) 2 ^A-*“ 2A 


-| 1/2 


or 


Ut.e. 

f/oo 


( U e T E \ ^ o 1 

^ ^ ' j + tan 2 A cos A (43) 


Design streamline boundary-layer method . 
One further method of evaluating the theoretical 
wake drag of a swept untapered airfoil will be consid- 
ered at this point. In the course of designing the con- 
toured liner (ref. 15) for the test of the suction airfoil, 
extensive theoretical airfoil boundary- layer calcula- 
tions (heretofore unpublished but used in the deter- 
mination of the liner contour) were carried out along 
the theoretical streamlines over the airfoil surface and 
into the wake. (See fig. 9(a).) Assuming that, for the 
design case at least, the cross-flow effects are small 
( v p << Up ), equations (33) and (34) become 


®uu — ®upu p + terms O(0 VpUp , 0 VpVp , 6 Up ) 
Ovu = 0 VpU p + terms O(0 VpUp , 9 VpVp , 6 * p ) 

Then equation (37) becomes 



+ terms 0{6 VpUp , 6 VpVp , 6* p ) (44) 

where the subscript “W.R.” denotes evaluation at the 
wake-rake position, X/c = 1.094. 

For this unpublished theoretical airfoil boundary- 
layer data computed along the streamlines, what is 
available from these computations is 

where U e = J u% + v\ and 

b'OO V 

Ve — ^00 = Uoo sin A, Uoo = Uoo cos A 


11 


Therefore 




and so equation (44) becomes 


c d,W 



The unpublished data used for equation (45), and the 
result of this computation, are presented in table IV 
for a range of boundary- layer transition locations. To 
compare this theoretical wake drag coefficient with 
that obtained experimentally, the assumed transition 
values of x/c were matched with those experimentally 
measured on the wing centerline at chord Reynolds 
numbers from 10 to 20 million. (See table V.) The 
wake drag coefficient computed at a chord Reynolds 
number of 20 million (the design-point Reynolds 
number at which the streamline boundary-layer com- 
putations were done) was scaled to lower experimen- 
tal Reynolds numbers on the basis of the well-known 
relationship for the turbulent wake of a flat plate 
(ref. 18, ch. XXI). 

c d,W V&c = Constant (46) 

These results are tabulated in table VI. The compar- 
ison of these theoretical wake drag coefficients with 
a typical experimental variation of wake drag coef- 
ficient with Reynolds number on the slotted LFC 
model using the conventional Baals and Mourhess 
method is presented in figure 9(b). 

Evaluation of Drag Coefficient Analyses 

Suction Drag 

Equation (13) was used to evaluate the suction 
drag on both surfaces of the swept LFC airfoil (fig. 1) 
over a range of chord Reynolds numbers at a constant 
Mach number. The integral in equation (13) was 
evaluated by a summation of the theoretical suction 


flow coefficient Cq (fig. 3 and tables 1(b) and (c)) 
over the airfoil for individual ducts and correspond- 
ing theoretical suction-pressure and temperature ra- 
tios to bring the exhaust suction air of the individual 
ducts back to free-stream velocity and pressure. Ta- 
ble I gives representative upper and lower surface val- 
ues of the design suction coefficients for several chord 
Reynolds numbers and design values of pressure and 
temperature ratios for individual suction ducts of the 
slotted airfoil. Also shown in table I are the design 
values of the calculated coefficient of suction drag, 
Cfi 5 , for upper and lower surfaces used to obtain the 
design suction drag. Equation (13) was used to ob- 
tain the suction drag coefficient for all the LFC tests 
in the 8-ft TPT. 

Wake Drag 

Five approaches to the computation of the wake 
drag on a swept 2-D airfoil have been presented. 
These are briefly summarized as follows: 

1. The conventional Baals and Mourhess method 
(ref. 2) for the calculation of the wake drag of a 2-D 
airfoil from static and total pressure measurements 
on a wake rake has been used to compute the wake 
drag for all of the LFC tests in the 8-ft TPT. It is not 
strictly correct for a swept airfoil. Estimates based 
on a brief tuft test indicate that sweep-induced cross 
flow at the wake rake is as much as 30°, and thus the 
wake drag from this method could be 5 to 10 percent 
high because of errors in total-pressure data due to 
misalignment between the local wake flow direction 
and the individual rake probes. 

2. The leading-edge-perpendicular method re- 
sulting in equation (29) is based on theoretical 
laminar-boundary-layer parameters computed from 
the component of free-stream flow that is perpendic- 
ular to the airfoil leading edge. These boundary-layer 
parameters could not be measured experimentally 
with the available instrumentation. The analysis of 
equation (29) indicates that the cross-flow terms are 
too large to ignore. 

3. The transformation to a coordinate system 
aligned with the free-stream direction results in 
equation (37). The cross- flow terms are less than 
0.2 percent of the total drag. The significant term 
in equation (37), the term involving the streamwise 
momentum thickness 0 UpUp , could be measured ex- 
perimentally but no such data were obtained during 
these tests. 

4. The semiempirical approximation of equa- 
tion (38), using the same experimental wake-rake 
data as the conventional Baals and Mourhess method, 
does not account for flow misalignment error that 
would cause the Baals and Mourhess drag to be 


12 



high. The two methods compare quite well at the 
lower Reynolds numbers, but begin to diverge for 
Reynolds numbers greater than 11 million, princi- 
pally in the lower surface contribution, probably be- 
cause the lower surface boundary layer begins to sep- 
arate at 85 percent chord at this Reynolds number 
(see ref. 8) and thus the shape factor H in equa- 
tion (38) becomes invalid. 

5. The design streamline boundary-layer wake 
drag method (eqs. (45) and (46)), which is a sim- 
plification of equation (37) by neglecting the cross- 
flow terms, yields values for wake drag variation with 
chord Reynolds number (fig. 9(b)) that compare well 
with the experimental values as calculated using the 
Baals and Mourhess method and thus tend to con- 
firm the liner design process. 

The data reduction program for the LFC tests 
in the 8-ft TPT uses the conventional Baals and 
Mourhess method for wake drag. All the error in 
this method results from misalignment between the 
wake-rake elements and the local wake flow. This 
error could be reduced if the rake could be placed 
further downstream of the airfoil trailing edge (but 
other design considerations made that impossible in 
this series of tests) or if the rake pitch could be varied 
at each test condition so that the true total pressure 
deficit was obtained on each probe. 

Tables II and III show that at the lowest Reynolds 
number for the theoretical boundary-layer compu- 
tations, R c = 8 x 10 6 , the three wake drag meth- 
ods (eqs. (29), (37), and (38)) using the theoretical 
boundary- layer characteristics yield, for full-chord 
laminar flow, values of = 2.5 x 10 -4 , 2.6 x 10~ 4 , 
and 3.0 x 10 -4 , respectively. These methods based 
on theoretical boundary-layer parameters agree well 
with each other, as did the methods based on ex- 
perimental wake data. Note that equation (38) 
was used both ways. However, the drag coefficients 
based on the theoretical full-chord laminar parame- 
ters are low compared with the experimental (Baals 
and Mourhess method) wake-rake drag coefficient of 


about 8 x 10“ 4 measured at R c = 10 x 10 6 even 
though those data were obtained with laminar flow 
to at least X/c — 0.95. It should be noted, however, 
that the experimental determination of 95-percent- 
chord laminar flow on the upper surface is based on 
a surface-mounted turbulence detector that was not 
in line with the wake rake but somewhat below it, 
as shown in figure 5. Extensive tuft studies indi- 
cate that even when this last downstream detector 
location is laminar, the wing upper surface flow at 
the trailing edge is very nearly parallel to the trail- 
ing edge and in the downward direction, so that the 
flow at the wake rake is probably contaminated to 
some extent by flow from nonlaminar regions fur- 
ther up on the wing. In addition, the patterns of 
laminar flow shown in reference 8 indicate that the 
lower surface flow immediately upstream of the wake 
rake was laminar only to about 60 percent (the flow 
above the wake rake was laminar to the end of suc- 
tion, 84 percent). 

Concluding Remarks 

The Langley Research Center has designed a 
swept, supercritical airfoil incorporating laminar-flow 
control (LFC) for testing at transonic speeds. Ana- 
lytical expressions have been developed and an eval- 
uation has been made of the drag, including suction 
and wake drag, by using theoretical design informa- 
tion and experimental data. 

The analysis shows that, although the sweep- 
induced boundary- layer cross-flow influence on the 
wake drag is too large to be ignored and there is not 
a practical method for evaluating these cross-flow ef- 
fects experimentally, the conventional unswept Baals 
and Mourhess method for wake drag computation 
used in the reduction of the experimental data yields 
wake drag that is at worst 10 percent too high. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 10, 1989 
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Appendix 

Evaluation of Total Suction Drag Coeffi- 
cient From Flow Quantity Measurements in 
Calibrated Suction Duct Nozzles 

The suction flow coefficient Cq of the individual 
suction ducts or chambers on the swept LFC airfoil 
may be defined as 


Cq,sc — 


Poo u oo^N c N 


(Al) 


where m sc j is the mass flow due to suction of the 
ith duct, which is located between successive chord- 
wise stations Sj and Si+i and has an average span bjy 
and chord cjq (measured parallel and perpendicular 
respectively to the leading edge). The duct dimen- 
sions include in every case half of the bulkhead wall 
thickness so that the sum of the duct areas is exactly 
equal to the projected airfoil surface area. Therefore, 
the total suction flow coefficient of the airfoil is equal 
to 


i=K 

(Cq,sc) total = 

i = 1 


=/' 


~~{p w )ws 
Poo^oo . 



(A2) 


where K is the number of suction chambers. Val- 
ues of Cq sc can be determined from flow quantity 
measurements of individual suction nozzles located 
in each duct (figs. 2(c) and 10). The following equa- 
tion for Cq sc is based on specially designed and cal- 
ibrated nozzles as previously used by Pfenninger and 
Groth in low-drag suction airfoil flight and wind- 
tunnel tests (ref. 19): 


Cq,sc ~ 


a N 




Esc^m 
Poo Qoo 


b N c N cos A yjl - (//F) 2 


(A3) 


where A ppr = Psc ~ PN 1S a required measurement. 
The term cos A in the denominator is required be- 
cause Cq sc is defined in equation (Al) as nondimen- 
sionalized by u 0 0 = Uqq cos A, but the q oo that is most 
convenient to use in equation (A3) is the free-stream 
dynamic pressure q oo = O-bpooU^. 

A sketch of the axisymmetric suction nozzle pre- 
viously used by Pfenninger and Groth (ref. 19) and 
considered herein is shown in figure 10. The nozzle 
flow coefficient apj is a function of Reynolds num- 
ber based on nozzle diameter, assuming geometri- 
cally similar nozzles, as determined by calibration. 


An empirical equation for the diameter of the nozzle 
from the beginning of the contraction to the throat 
is (for 0 < x < 2 dpj) 


d(x) = dpj 





( A4 ) 

where dp? is the nozzle throat diameter and the 
exponent n = 2.545. This equation fits the nozzle 
coordinates from the entrance plane ( x = 0) to the 
throat where ppi is measured (x — 2 djy), to within 
0.01 in. for d^ = 1 in. In the 8-ft TPT LFC 
experiments, 12 sizes of this nozzle were used, with 
dps = 0.188 to 0.875 in. 


The data reduction procedure obtains the value 
of c*jv from an empirical equation fitted to the nozzle 
calibration data. The fit of the empirical equation to 
the data is within the error band of the calibration 
data. The empirical equation used in the LFC data 
reduction is 


a N (R d ) = 0.927 - 0.0339787 


x 


v /l + 4(log 10j R d - 4)2-1 


+ 0.075(log 10 R d ~ 4) (A5) 


where R d is the Reynolds number of the nozzle 
throat flow based on the nozzle throat diameter 
dp}. This equation was fitted to data in the range 

10 3 <R d < 10 5 . 

The suction drag may finally be expressed (eq. 13)) 
as follows: 


cos A 
0.2MI 

CM 

0.2MI 

+ T l: sc [ 




/ \ 
(Poo \ 1 

1 

) d(-) 

Too . 

\Psc) 

- 

i ^ c 


The experimental accuracy to which the suction 
drag coefficient can be determined for a single duct 
is primarily dependent upon the accuracy to which 
A p (the drop in pressure from the duct to the noz- 
zle throat) can be measured. Inspection of equa- 
tion (A3) for the suction flow coefficient reveals that 
Cq,sc is proportional to (A Pat) 0 ' 5 with all other quan- 
tities constant, except for the effect of which 
is a function of nozzle throat Reynolds number and 
is determined from calibration. Therefore, in low- 
drag experiments involving determination of the suc- 
tion drag as part of the total drag, when equation 
(A6) is applied, the data accuracy required should be 
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analyzed and a suitable measurement technique (in- perpendicular to the wing leading edge, while the de- 

strumentation) selected. Note that the sweep A ap- sired suction drag coefficient c^ s is in the free-stream 

pears in equation (A6) (and eq. (13)) because Cq, sc is direction, 

defined in terms of the velocity Uqq = Uqq cos A, 
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Table I. Representative Theoretical Values of the Suction 
Drag Coefficient and Suction-Chamber Parameters on 
the Upper and Lower Slotted LFC Airfoil Surfaces 

[Moo = 0.82, A = 23°] 


(a) Suction coefficient and suction drag coefficient 
as integrated from table 1(b) parameters 


R c 

Cq 

Cd,s (eq- (13)) 

Upper surface 

8 x 10 6 

3.848 x 10~ 4 

5.164 x 10~ 4 

20 

2.970 

4.051 

40 

2.464 

3.408 

Lower surface 

8 x 10 6 

6.568 x 10~ 4 

7.054 x 10~ 4 

20 

5.142 

5.492 

40 

4.319 

4.589 




Table I. Continued 


[R c = 20 x 10 6 


(b) Upper surface suction-chamber pressure and temperature ratios 
and suction-coefficient variation with X/c 



[ x 10“° 
! x icr 3 
l x 10- 2 
I x 10 -2 
! x 10~ 2 
I x 10~ 2 
I x lO -1 
x 10 _1 
I x KT 1 
! x 10 _1 
l x 10 -1 
; x to -1 

’ x 10 _1 
I x KT 1 
x icr 1 
’ x KT 1 
’ x KT 1 


0.6113 

.6113 

.6345 

.6345 

.6345 

.6345 

.6445 

.6445 

.6445 

.6498 

.6498 

.6530 

.6530 

.6530 

.6530 

.6577 

.6577 

.6607 

.6607 

.6653 

.6653 

.6653 

.6653 

.6695 

.6695 

.6695 

.6722 

.6722 

.6722 

.6722 

.6727 

.6727 

.6727 


0.9570 

.9570 

.9586 

.9586 

.9586 

.9586 

.9596 

.9596 

.9596 

.9603 

.9603 

.9608 

.9608 

.9608 

.9608 

.9613 

.9613 

.9618 

.9618 

.9624 

.9624 

.9624 

.9624 

.9630 

.9630 

.9630 

.9634 

.9634 

.9634 

.9634 

.9638 

.9638 

.9638 


-0.1200 x 
-.1160 x 
-.1130 x 
-.1080 x 
-.1040 x 
-.1060 x 
-.1040 x 
-.1020 x 
-.1000 x 
-.9900 x 
-.9800 x 
-.9800 x 
-.9400 x 
-.9300 x 
-.9300 x 
-.9200 x 
-.9000 x 
-.8800 x 
-.8600 x 
-.8200 x 
-.8100 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8000 x 
-.8200 x 








Table I. Continued 


(b) Concluded 
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Table I. Continued 


[Re - 20 X 10 6 ] 


(c) Lower surface suction-chamber pressure and temperature ratios 
and suction-coefficient variation with X/c 



0.1000 x 
.6228 x 
.2102 x 
.4595 x 
.8159 x 
.1277 x 
.1842 x 
.2504 x 
.3204 x 
.4122 x 
.5077 x 
.6128 x 
.7280 x 
.7925 x 
.8725 x 
.9525 x 
.1020 
.1045 
.1065 
.1080 
.1095 
.1115 
.1145 
.1230 
.1291 
.1304 
.1420 
.1490 
.1515 
.1535 
.1550 
.1565 
.1580 
.1600 
.1657 
.1728 
.1849 
.1969 
.2093 
.2224 
.2361 
.2503 


0.1143 x 10 1 

0.9793 

-0.1000 x 10 -3 

.1153 x 10 1 

.9804 

-.3000 x 10 -3 

.1153 x 10 1 

.9804 

-.4000 x 10~ 3 

.1153 x 10 1 

.9804 

-.6000 x 10~ 3 

.1103 x 10 1 

.9843 

-.9500 x 10~ 3 

.1103 x 10 1 

.9843 

-.2500 x 10~ 2 

.1103 x 10 1 

.9843 

-.4700 x 10 -2 

.1031 x 10 1 

.9884 

-.3600 x 10 -2 


.1114 x 10 1 
.1114 x 10 1 
.1114 x 10 1 
.1114 x 10 1 


.1900 x 10" 2 
.3400 x 10 -2 
.3500 x 10~ 2 
.3300 x 10“ 2 


.1500 x 10" 2 
.1500 x 10" 2 
.1000 x 10~ 2 


.7200 x 10~ 4 
.7200 x 10" 4 
.7200 x 10~ 4 


I 



Table I. Continued 


(c) Continued 
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Table II. Wake Drag Computations Prom Representative Theoretical Trailing-Edge 
Boundary-Layer Values on the Upper and Lower Surfaces of the Slotted LFC Airfoil 


[Moo = 0.82, A = 23° ] 

(a) Theoretical trailing-edge boundary-layer parameters with full-chord laminar flow 


R c 

“e,T.E. 

Woo 

U e , T.E. * 

®uu 

c 

$vu 

C 

61 

£ 

1 

c 

Upper surface 

1 

K 

0.9168 

0.9300 

1.3820 x 10~ 4 
.8740 
.6180 

0.9484 x 10~ 4 
.5998 
.4241 

4.072 x 10“ 4 
2.575 
1.821 

2.782 x 10~ 4 

1.760 

1.244 

Lower surface 

8 x 10 6 
20 
40 

0.8293 

0.8576 

0.8662 x 10~ 4 
.5479 
.3847 

0.8504 x 10~ 4 
.5378 
.3803 

1.842 x 10" 4 
1.165 
.824 

2.203 x 10~ 4 
1.393 
.985 


t Equation (43). 


(b) Representative wake drag coefficients computed with equations (29) and (38) 


Rc 

c d,W 

(eq. (29)) 

2nd term,J 
percent 
(eq. (29)) 

c d, total 

(eqs. (13) and (29)) 

c d,W 

(eq. (38)) 

c d, total 

(eqs. (13) and (38)) 

Upper surface 

8 x 10 6 

1.771 x 10" 4 

14 

6.935 x 10~ 4 

2.072 x 10 -4 

7.236 x 10~ 4 

20 

1.120 

14 

5.171 

1.310 

5.361 

40 

.792 

14 

4.200 

.926 

4.334 

Lower surface 

8 x 10 6 

0.892 x 10“ 4 

22 

7.946 x 10~ 4 

1.002 x 10" 4 

8.056 x 10~ 4 

20 

.564 

22 

6.056 

.634 

6.126 

40 

.399 

22 

4.988 

.448 

5.037 


f The 2nd term consists of multiples of 0 VU . 


23 

























Table III. Wake Drag Computations From Representative Theoretical Trailing-Edge 
Boundary-Layer Values on the Upper and Lower Surfaces of the Slotted 
LFC Airfoil Based on Coordinate Transformation of Equation (32) 


[Moo = 0.82, A = 23°] 


(a) Theoretical trailing-edge boundary-layer parameters with full-chord laminar flow 


Rc 

@UpUp 

c 

@vpUp 

c 

@VpVp 

c 

5 

c 

C 

Upper surface 

8 x 10 6 
20 
40 

1.3350 x 10" 4 
.8444 
.5971 

3.230 x 10 -6 

2.040 

1.440 

3.061 x 10~ 5 

1.936 

1.369 

3.836 x 10~ 4 

2.426 

1.716 

4.979 x 10“ 5 

3.149 

2.227 

Lower surface 

8 x 10 6 
20 
40 

0.8864 x 10~ 4 
.5606 
.3964 

-0.690 x 10" 6 

-.440 

-.310 

0.745 x 10 -5 
.471 
.333 

1.992 x 10 -4 
1.259 
.891 

-1.45 x 10" 5 

-.915 

-.647 


(b) Representative wake drag coefficients computed with equation (37) 


R c 

c<t,W 

(eq- (37)) 

3rd term,t 
percent 
(eq. (37)) 

c d, total 

(eqs. (13) and (37)) 

Upper surface 

8 x 10 6 

1.553 x 10~ 4 

0.2 

6.717 x 10 -4 

20 

.982 

.2 

5.033 

40 

.694 

.2 

4.102 

Lower surface 

8 x 10 6 

0.985 x 10 -4 

0.08 

8.039 x 10 -4 

20 

.623 

.08 

6.115 

40 

.441 

.08 

5.030 


t The 3rd term consists of multiples of 6y p . 
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Table IV. Design Streamline Boundary- Layer Profile 
Functions Used to Compute Wake Drag From Equation (45) 


Xtr 

ZW.R. 

Ue/Uoo 

6* 

Up 

@u p u p 

c d,W 

Upper surface 

0.90 

1.094 

0.9179 

4.396 x 10“ 4 

3.045 x 10~ 4 

4.183 x 10~ 4 

.80 



8.062 

5.433 

7.451 

.76 



8.752 

5.870 

8.048 

.74 



9.517 

6.352 

8.706 

.67 



13.58 

8.845 

12.10 

.60 



16.31 

10.48 

14.33 

.55 



20.11 

12.68 

17.32 

.50 



22.09 

13.80 

18.83 

.45 



25.07 

15.46 

21.08 

Lower surface 

0.84 

1.094 

0.9191 

2.703 x 10~ 4 

1.962 x 10“ 4 

2.715 x 10“ 4 

.75 



3.265 

2.360 

3.265 

.65 



4.821 

3.477 

4.809 

.55 



8.502 

6.089 

8.419 

.44 



13.49 

9.529 

13.16 

.20 



26.91 

18.39 

25.36 

.17 



27.84 

18.99 

26.18 

.07 



28.78 

19.62 

27.05 

.05 



28.78 

19.62 

27.05 


Table V. Variation of Boundary-Layer Transition Location Along Streamline 
Through Wake Rake and Wake Drag, With Reynolds Number 
for Optimum Slotted LFC Suction Data 


Rc 

Upper 

Lower 

Total 

Forward 

X/c 

Aft 

X/c 

Forward 

X/c 

Aft 

X/c 

c dW 

10 X 10 6 

> 0.90* 

1.00 

0.84 

0.93 

6. 7-8.3 x 10~ 4 

11 

> .90* 

1.00 

.85 

.97 

7.0-8.2 

12 

.74 

1.00 

.89 

.95 

10.2-11.3 

13 

.60 

.76 

.64 

.75 

14.4-15.3 

14 

.53 

.67 

.54 

.74 

20.1-20.4 

15 

.55 

.67 

.17 

.44 

38.8-39.2 

16 

.55 

.67 

.17 

.44 

41.8 

18 

.52 

.67 

.17 

.44 

42.4-42.7 

20 

.52 

.67 

.04 

.44 

41.8-43.4 


t Default transition at flap hinge. 




Table VI. Theoretical Wake Drag Variation With R c Based on Streamline Boundary- Layer 
Analysis at R c = 20 x 10 6 Rescaled to Test R c on the Basis of c^w Rc^ = Constant 
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t Default transition at flap hinge. 

* Default transition at end of suction. 




Tollmien-Schhchting 


Supersonic 






: LFC airfoil laminar i 


1 region). 
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Contoured top 
tunnel wall 


Flap hinge, 
X/c = 0.89 


Turbulent wedge V 'T 
X / / / //Laminar/ / / /xt7 

w|p| 

* /" Turbid 


End suction 
X/c = 0.842 


(b) Lower-surface test region. 


Figure 2. Continued. 




Slotted skin-plenum-duct design 



Porous skin-flute-duct design 



(c) Suction system design cross sections. 


Figure 2. Concluded. 















TOP VIEW 


Hinge Trailing Total tube 
line edge face 




SIDE VIEW 



Figure 5. Schematic showing orientation of wake rake to LFC airfoil in the 8-ft TPT. 


Pt 
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Pt 

P t , 00 


1.12 
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.90 

.88 

.86 




O Rc = 10 x 10 6, p t = 751.3 psf 
□ R c = 1 1 x10 6,p’ =824.0 psf 
OR c =12x106, p J’ co = 900.9 psf 


3 4 

y, in. 



O R c = 1 3 x 1 0 6, p. = 975.0 psf 
□ R c =14x10 6 , p = 1 050.9 psf 
OR c = 15 x 10 6, pj’”= 1125.8 psf 


y, in. 

Figure 6. Typical variation of wake-rake static and total pressure ratios across wake (y = 0 on lower 
surface side) of swept LFC airfoil at design Mach number. 


1.00 


O Rc= TO X TO®. Pt co = 7513 P sf 
□ R c = 1 1 x10 6,p' =824.0 psf 
OR c = 12x106, p ' = 900.9 psf 


O Rc= 13 x 10 6, p. = 975.0 psf 
□ R c =14x10 6, p!’ = 1050.9 psf 
OR c = 15 x 10 6, p{’”= 1125.8 psf 


y, in. 

? igure 7. Typical variation of wake velocity ratio across wake (y = 0 on lower surface side) of swept LFC 
airfoil at design Mach number. 



Rc 

Figure 8. Comparison of typical wake drag coefficient variation with chord Reynolds number as computed 
from equation (38) with that computed by the Baals and Mourhess (ref. 2) 2-D method. 
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c d,W 25 


/ °/ 
oy 


O Experiment (slotted U.S.) 
Baals & Mourhess wake drag 

Design theory eq. (45) 

transition forward 

Design theory eq. (45) 

transition aft 


10 11 12 13 14 15 16 17 18 19 20 21 X10 6 

Rc 

(b) Comparison of theoretical wake drag with experiment. 

Figure 9. Concluded. 
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